Bistability and Hysteresis in the Sliding Friction of a Dimer 
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The sliding friction of a dimer moving over a periodic substrate and subjected to an 
external force is studied in the steady state for arbitrary temperatures within a one- 
dimensional model. Nonlinear phenomena that emerge include dynamic bistability 
and hysteresis, and can be related to earlier observations for extended systems such 
as the Frenkel-Kontorova model. Several observed features can be satisfactorily ex- 
plained in terms of the resonance of a driven-damped nonlinear oscillator. Increasing 
temperature tends to lower the resonant peak and wash out the hysteresis. 

PACS numbers: 81.40.Pq,46.55.+d 



1. INTRODUCTION 

The friction experienced by atoms, small molecules, and adlayers moving over substrates 
is an active topic of current research pj. One reason is the desire to understand, at a 
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fundamental level, the origin of friction. The other reason is the wish to acquire expertise 
in developing nano-devices like nano-motors, nano-wires, and nano-probes. Besides, sliding 



'aces, such as diffusion of atoms and 
, 10] over periodic substrates. That 



friction is related to other atomistic processes at sur 
molecules P, [], U, 0] , and motion of long chains P, Q, E , 
microscopic sliding friction can exhibit nonlinear behavior depending on the sliding regime 
has been shown in the recent theoretical literature. Strunz and Elmer jfj] have studied in 
detail the nonlinear friction of the Frenkel-Kontorova model and identified the origins of 
such friction as being the resonance of the sliding velocity with the internal vibration modes 
of the chain, and the formation of kinks. The latter phenomenon was also studied by Braun 
et a, m. More r ecen tly , FuS co and Fasoiino 0Q nave identified tfie SaI ne resonance 
phenomenon in the friction of a smaller object, a dimer moving over a periodic substrate. 
Gongalves et al. [3] have analyzed a closely related system in the relaxation regime, i.e. in 
the absence of external forces, and have given a simple physical explanation of the observed 
results. Persson |lJ| has predicted that a friction force proportional to v ~ 3 is to be expected 
(in addition to the linear one) in the large-force regime for a sliding system of any size, 
ranging from a single atom to an infinite chain. In the present paper, we restrict our study 
to the steady state friction of a dimer, but investigate all force regimes. We identify several 
separate regimes exhibiting resonance, bistability, and hysteresis. We provide a simple 
explanation for these phenomena in terms of the resonance of a driven-damped oscillator, 
and comment on how they depend on the substrate corrugation, damping and temperature. 

2. MODEL AND SIMULATION RESULTS 

At zero temperature, the equations of motion for the two particles constituting the dimer 
sliding over a periodic potential, in the presence of external force F, are Q| : 

• ,/ x 27TM . f 27TXi\ 

mx\ + m^fXi — k[x 2 — X\ — a) = — — sm I — - — + t 



2tTU ( 1llX2\ „ 

mx 2 + m^Xi + k\x 2 — X\ — a) = — — sm — - — + F, (1) 

b \ b J 

where x\^ are the coordinates of the two particles each of mass m, and k, a, b, u are, 
respectively, the spring constant, equilibrium length of the dimer, wavelength of the substrate 
potential and half the amplitude of the potential. We integrate these coupled equations 
numerically using the algorithm of Verlet modified to allow for velocity-dependent forces 
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FIG. 1: Force- velocity (f /"f versus v, f = F/m) relation obtained from numerical integration in the 
steady state. Both f/'y and v are expressed in units of bu>o. Parameter values are u = 0.038fc6 2 , 7 = 
(2/3)a;o, and a/b = 1/2. Notice that the maximum friction at resonance occurs when 2irv/b uo, 
therefore for abscissa values which are around (27r) _1 . 

The underlying characteristic physical quantities in this system are the equilibrium dimer 
length a, the substrate wavelength b, the free dimer characteristic time l/u>o — a/ m/2k, and 
an energy that describes the dimer oscillation such as kb 2 . In our numerical integrations we 
use a time step At equal to 0.03UQ 1 . 

The procedure is as follows: for a fixed value of F we obtain the center of mass velocity 
v, averaged over several thousand time steps in the steady state. Repeating the procedure 
for several hundreds of different values of the force, we construct a characteristic curve of 
force versus velocity, F/{m^) versus v, (see FigJTJ). 

The following features are evident from FigQ]as emerging from our numerical simulations: 

1. Existence of a static threshold: there is a minimum value of the externally applied force 
required to make the dimer slide. It arises from the fact that, at zero temperature, 
the substrate potential has to be overcome. 

2. Linear behavior in the asymptotic large-force regime: for sufficiently large forces, the 
dimer slides at velocities high enough to make the substrate potential a negligible 
perturbation. The linear asymptote is shown in Figdas the dashed line. 
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3. Nonlinearity in the friction exhibiting a maximum: this arises from resonance effects 
when the washboard frequency 27tv/b is close in value to the dimer frequency ujq = 
^2k/m. 

The first two of these features are expected in light of previous reports in the literature. The 
third appears to be interesting. Previous related reports have been on the driven diffusion 
of a dimer [l^j], and of a Frenkel-Kontorova model in the presence of an external force 
A thorough analysis is however lacking. We attempt to provide such an analysis below in 
the specific case of the dimer. We relate our results to observations made by other authors 
in the context of the Frenkel-Kontorova model, and suggest that the essential features of 
the nonlinear friction of the infinite linear chain can be understood in terms of the dimer 
dynamics we describe. 



3. SIMPLE THEORETICAL CONSIDERATIONS 



A transformation of the coordinates of the two dimer masses to the center of mass co- 
ordinate x + = (x 2 + xi)/2, and the internal coordinate x_ = (x 2 — £i)/2, converts Eqs.JT]) 
to 

2nu . {2nx + \ {2nx_\ . 
x + + 1 x + = —s^—jcos^— j+F/m 

k 2nu ( 2ixx\ f 2nx + \ 

a:_ + H (2a;_ — a) = — — sm — - — cos — ; — . (2) 

m mb \ b J \ b J 

If in the last equation we define £ = — 1 = X2 ~ Xl — 1, we get 

as in the analysis of Ref. Contrary to that analysis, however, here our interest lies in 
the steady state reached after the application of the external force. In that state, the center- 
of-mass velocity v + (t) oscillates around a constant value v. Considering only situations in 
which Av(t) = \v+(t) — v\ < v, let us neglect Av(t) and decouple the equations. The internal 
coordinate then satisfies 

J + + * = sm (t (1 + «) C08(w ' ) - (4) 

which describes a damped nonlinearly driven oscillator. The natural frequency is ujq — 



^J2k/m, the damping is 7, and the driving frequency uj is proportional to the constant 
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component of the center of mass velocity: uj = 2nv/b which is the so-called washboard 
frequency. 

3.1. Linear analysis in zeroth order 

Eq.(j3I) cannot be solved analytically because of the nonlinearity in the sine term. The 
simplest approximation, valid in zeroth order in powers of £, leads to the equation of a 
driven-damped linear oscillator: 

^ + 7 ^ + ^ = ^ sm lTj cos( " t) ' (5) 

with an exact solution in the steady state, 



Amu . /7ra\ 1 , 

£(*) = — - sin — = cos (ut - 5) , (6) 



where 5 is the phase angle given by tan(<5) = ^rzji- In order to compute the dependence 
of the friction on the center-of-mass velocity, we must resort to the power balance condition 
which, in terms of the averaged center-of-mass and internal velocity, can be written as 

F < v + >= rwy < v\ > +mj < v 2 _ > . (7) 

Defining / = F/m, and using the previously made assumption that < v\ >~ v 2 , leads to 

<v 2 _> 

J=iv + 1 . (8) 

v 

We see that, generally, the steady-state friction the center of mass of the dimer experiences 
is nonlinear in the velocity. Using the steady state solution for £(t) as given by ©, we 
calculate 

i /n \ 2 2 

2 \l 2ttu \ . 2 (-na\ UJ 



< v l >= - sin 2 ^ — 9 . (9) 

Substitution of < v 2 _ > given by Eq.fjnj in Eq.JHJj yields 

f --v = \ ( sin 2 C4) V - . (10) 

7 2 \m) V b ) ( w 2 -(^o/2vr) 2 ) 2 + t; 2 (&7/27r) 2 V ' 

The right hand side of Eq. lfTUl) focuses on the nonlinear component of the friction. In FigJ21 
we plot Eq. (fT0|l together with the exact results from the simulations. We see that our zeroth 
order theory captures the essence of the resonance behavior, approximates well the location 
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FIG. 2: Resonance effect in the nonlinear friction of the dimer and comparison of the simulations 
with the zeroth order theory. Plotted is the nonlinear part, f /j — v of the characteristic curve for a 
dimer sliding over a periodic substrate. Circles are from numerical integration at the steady state. 
The solid line is the zeroth order theory for the uncoupled equation for £. f/j and v are both in 
units of biOQ. The parameters in this case are u = 0.028/c& 2 ,7 = (2/3)uo, and a/b = 1/2. 

of the peak, is excellent quantitatively for larger velocities, and fails only to describe the 
low-velocity threshold. 

In the limit of high velocities (v >> buo / 2tt , b'-f / 2tt) , the leading term in the resonance 
denominator is quartic in v, resulting in the following nonlinear friction: 



This agrees with the 1/v 3 result of the sliding friction of the purely internally damped dimer 
studied recently In the other limit, in the low- velocity regime (v « bu>o/2ir), the 

denominator of the right hand side of Eq.([T0]) is independent of v (it is of fourth order in 
c<j ), so f /j — v in this case becomes linear: 



Here v s = bcu /2ti is the sliding velocity near resonance. Thus, in the limit of typical 
velocities of experiments, like the microbalance experiment, we recover the linear regime 




(11) 




(12) 




represents the part of the friction directly contributed by the 
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interaction with the substrate potential. However, the threshold effect in the dimer problem 
at T = prevents this regime from being observed. For an extended object and/or at 
T > 0, where the threshold could be vanishing, it would be possible in principle to observe 
this regime. 

3.2. Parametric oscillator in the first order 

The zeroth order description fails when a/b is commensurate, i.e. when takes integer n 
values. In this case the zeroth order approximation cannot be applied and one is forced to 
go to the next order, since neglecting £ in the sine terms predicts an erroneous (vanishing) 
nonlinear friction. In such cases, and if £ << 1, 

sin + 0) = (-1)" sin(nvrO » (-l) n rnc£. (13) 

Substitution converts the driven-damped harmonic oscillator equation (J3| into the equation 
for a parametric oscillator 

w +7 i + H + ^sr cosM) ) e = - (14) 

For this equation an exponential increase of £ is expected in an instability window around 
uj = 2u . Thus, in this regime, < vt. > would increase indefinitely and the friction force 
would be infinite, which is of course unphysical. In fact, in the full system Eq.J2}, the 
coupling between the center-of-mass and internal motion drives the center of mass out of 
the instability window characterizing the parametric oscillator, and this is enough to make 
the increase of £ saturate |ll|, Q], yielding a finite value of < vt >■ Besides, if £ increased 
because of the resonance, the assumption (|T3|) that leads to the parametric oscillator equation 
AH} would not be valid. 

The zeroth order approach predicts that the shape of the resonance is proportional to 
sin 2 (7ra/6). Accordingly, there are two extreme cases worth discussing: a/b = 0.5 + nn 
in which the resonance is maximum, because the counter phase movement of the particles 
makes this the most energy effective case; in contrast a/b = nn is the less effective one, since 
the in phase movement does not excite the internal vibration, thus yielding purely linear 
friction. The latter case is precisely the one discussed in this section. Figure 03 illustrates the 
previous discussion showing how the resonance is actually affected by the commensuration 
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FIG. 3: Force-velocity relation obtained from numerical simulation for different values of a/b. 
Notice how the resonance is attenuated as a/b takes on increasing values from 0.5 to 1 and how the 
case a = b is very close to the linear regime. The difference can only be appreciated in the inset 
where f/y — v versus v is plotted on log-log scale, f/j and v are both in units of bu>o. Parameter 
values are u = 0.0056A:6 2 ,7 = (l/6)w - 

ratio a/b. Those are results of the numerical integration of Eqs.JH, where we can see that 
the zeroth order approximation is very good indeed. In the case a/b = 1, although the 
friction is not exactly the linear one, there is no resonance at all and the f-v characteristic 
is very close to the linear one. In the inset however, we can see that there is nonlinear friction 
which goes asymptotically v~ 3 to the linear regime, due to the parametric resonance. 

3.3. A non-standard description in higher order 

A non-standard approximation procedure may be developed through an iterative pro- 
cedure by noticing that, as the value of an arbitrary variable y increases, the expression 
cos(/« cos(y)), where n is a constant, oscillates around the value Jo( K ), whereas sin(/« cos(y)) 
oscillates around 0, Jo being the Bessel function of order 0. We have seen that, considered 
in zeroth order, Eq. predicts that £(t) oscillates sinusoidally with frequency to (see ©). 
Writing p = 7ra/6, we may use the fact that 

< sin (p(l + £(t))) >^< sin(p) cos(p£(£)) >= sin(p) J (pA sm(p)Z(uj)) , (15) 
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where A = and Z{oS) = —, =5 . This procedure can be followed iteratively to 

whatever degree is desired. The successive approximations to the steady state are thus 

f°(t) = Asin(p)Z(o;)cos(u;t-5) 

^(t) = Asm(p)Z(uj)J (pAsm(p)Z(uj))cos(ujt- 5) 

£ 2 (t) = Asm(p)Z(Lu)J (pAsm(p)Z(uj)J (pAsm(p)Z(Lu)))cos(ujt- 5)... (16) 

Generally, this may be expressed by defining rj = A sm(p)Z(uj) and writing the approxima- 
tion as 

f (t) « vJoipnJoipnMpnJoipn ...)))) cos(u;t - 5). (17) 

The spectrum predicted through this approximation is seen to be proportional to 
^ 2 V 2 Jo(P r }JoipvJoipvJoipV- ■•)))) to whatever degree of approximation one requires, except 
for p ps n7r, in which case the approximation is not valid, for the same reason explained in 
the previous section. 

The left panel of FigQ] shows a comparison of different orders of this non-standard ap- 
proximation with the numerical solution of Eq.(J3J), shown as a solid line. The zeroth order 
approximation in this procedure, shown as dotted line, overestimates the height of the reso- 
nance peak, while the first order approximation, shown as dashed-dotted line, underestimates 
it. In the second order, shown as dashed line, our procedure is already able to essentially 
coincide completely with the numerical solution. It is important to realize that our non- 
standard procedure does very well within the second order when viewed as an approximation 
to Eq.flU), for which it has been developed, rather than to the original Eqs.((H). The right 
panel shows the comparison of the numerical solution of Eq. (j3jl ; in which the center of mass 
velocity is a constant, with the simulation based on the original Eqs.JTJ). One can see a 
low-velocity departure arising from the static threshold and a high-velocity departure aris- 
ing from bistability. Such differences arise from the fact that in Eqs.fH the center-of-mass 
velocity is itself decided by the dynamics, therefore not being a free parameter as in Eq.(JH). 
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FIG. 4: Validity of a non-standard approximation procedure. Plotted is: f/j — v versus v, both in 
units of bujQ. (a) Comparison between the numerical simulation of Eq.Q for the internal coordinate 
£, with the nonstandard approximation in various orders, showing excellent convergence within the 
2 nd order, (b) Comparison between exact solution of system (QJ with the solution of Eq.(@J), in 
which the center of mass velocity is a free parameter; notice that except for the discontinuities the 
agreement is excellent. Parameter values are u = 0.016A:6 2 ,7 = (l/4)w , and a/b = 1/2. 

4. FURTHER NONLINEAR RESULTS 
4.1. Bistability and Hysteresis 



The nonlinearities present in our system give rise to bistability. This can be seen in 
FigJHl where it is clear that one value of the force can correspond to two distinct values of 
the velocity, in a certain region. In order to gain deeper insight into this issue, we plot in 
FigOHb) the prediction for the characteristic curve based on expression fTn| as u increases. 
We see features typical of bistable systems such as encountered in the pressure- volume (p- 
V) diagram of a van der Waals gas. The bistability of the dimer is directly related to 
hysteresis. FigEl shows the presence of two regions, identified by the dashed lines, where 
hysteresis occurs when the force is first increased from / = and then is decreased in small 
steps down to / = 0. The first hysteresis at low velocities is due to the bistability between 
the locked and the running state and has a static origin, being associated with the energy 
threshold that the particle has to overcome in order to move. The same kind of hysteresis 
is also found in the underdamped monomer On the other hand, the second hysteresis 
at intermediate values of v is of purely dynamical nature, and it is related to the bistability 
between two running states in that region. Let us discuss the latter hysteresis in more detail. 
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FIG. 5: Force- velocity curves f/j versus v (both in units of buj$) for different values of the corru- 
gation amplitude u expressed in units of kb 2 . Parameter values are: a/b = 0.5,7 = 0-/^)^0- (a) 
Simulation results (numerical integration of Eqs.JU). (b) Zeroth order approximation (Eq. (|lf)|l). 
Notice in both panels that the critical value of u (i.e. the value at which the bistability develops) 
is around 0.004/c6 2 for the chosen values of parameters. 

Five regions, denoted by roman numbers, can be identified there : I) the mechanically stable 
region where v increases with F up to a point where the force begins to be bivaluated, II) the 
metastable region where v increases with F up to the local maximum, III) the mechanically 
unstable region where v decreases with F up to the local minimum, IV) the metastable 
region, where v increases with F from the local minimum up to a point where the force ends 
to be bivaluated, and V) the mechanically stable region where v increases again with F and 
it is monovaluated. Region I) and V) are both obtained either increasing the force from 
zero or decreasing it from a high value. Region II, however is accessible while increasing 
the force from zero but not when the force is decreased from a high value. Complementary, 
region IV is obtained while decreasing the force from a high value but not when the force 
is increased from zero. Region III) is mechanically unstable, i.e the velocity decreases when 
the external force is increased. For a given value of the external force, the system goes to 
a stable point, therefore when the force goes beyond the value at the local maximum there 
is a jump from region II to region V. Conversely when the force is decreased from region V 
the system enters the metastable region IV until at the local minimum it jumps to region I. 
Notice that, except for the unstable region III, all the hysteresis features of the system are 
excellently reproduced by the solution of Eq.(|I|). 

FigEI shows the behavior of the center-of-mass velocity for two different, but close, val- 
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FIG. 6: Hysteresis in the force-velocity relation. Symbols are results from simulations of the system 
Eqs.Q: (□) increasing the external force, (<) decreasing the external force. The solid line is the 
result of numerical integration of the Eq.(@J). I and V indicate the stable regions (overlap of the 
two symbols), II and IV are metastable regions, while III denotes the unstable region. The arrows 
indicate the forward jumps (— >) and backward jumps («— ). Plotted is f/j versus v (expressed in 
units of bu>o). Parameter values are u = 0.023M> 2 ,7 = (l/6)w 0) and a/b = 1/2. 

ues of the external force near the region of bistability. In the early stages of the dynamics 
the velocity is practically the same in the two cases, but after some time the velocity cor- 
responding to the lower force attains a steady state value that is much smaller than the 
one corresponding to the larger force. This clearly illustrates the dynamical origin of the 
bistability. Furthermore, the oscillations of v + in the steady state are much smaller for the 
larger force. The bistable behavior of the dimer critically depends on the parameters u, u 
and 7: For fixed u and 7, it is observed when u exceeds a critical value, which can be 
estimated in the framework of the zeroth order approximation as the value for which the 
local maximum and the local minimum in the velocity-force characteristic, given by Eq. ljlOj) 
by adding the linear term in v, coincide. This is related, in the van der Waals analogy, to 
the region of coexisting phases. The interplay between linear and nonlinear friction due to 
resonance gives rise to the bistability, as does the interplay between attractive and repulsive 
terms in the van der Waals gas. In this sense, Eq. (fT0|) might be regarded as an equation of 
state of the system. Both equations are however only approximate descriptions of the real 
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FIG. 7: Bistable behavior exhibited in the time evolution of the center-of-mass velocity for slightly 
different values of the external force / = F/m. Velocity v + is in units of bu®, time t in units 
of 0Jq 1 , and / in units of kb. Notice that, up to time t = QOoJq^ , the two curves follow almost 
the same dynamics, but then separate, respectively, to a low and high limit. Parameters are 
u = 0.011M> 2 , 7 = (l/6)w , and a/b = 1/2. 

system and cannot predict what happens in the transition region, where a coexistence of 
two different states is found. 

Since the hysteresis is intimately linked to the bistability, it appears only for large u 
(and/or small 7), as discussed above. 



4.2. Effects of non-zero temperature 

It is natural to inquire into the effects of finite temperature in our system. We solve 
Eqs.(IH) by adding random forces representing the thermal interaction with the substrate, in 
the Langevin approach: 



, , N 27TM . 

mx\ + rwyxi + k[X2 — x\ — a) = —— sm 



27TM /27TX 2 \ „ _ 

mx2 + 7727x2 — k(X2 — x\ — a) = sm I — - — J + R 2 + F, (18) 
in which the stochastic forces Ri^ satisfy the conditions 

< R(t) > = 
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FIG. 8: Temperature effects on the force- velocity relation. Plotted is the f/*y versus v relation 
as temperature T (expressed in units of u/ks) is varied, f/'y and v are in units of bvo, and 
parameters are u = 0.011A;6 2 ,7 = (l/6)wo, and a/b = 1/2. (a) Curves for increasing force, the 
inset shows the details in the static threshold region. The bistability is seen to disappear as T 
increases smoothening the characteristic f-v curve, (b) Forward and backward curves to emphasize 
the hysteresis dependence on temperature: the area of the hysteresis loop decreases with T and 
disappears for temperatures close to fcgT = 0.005n. 



where fee is Boltzmann's constant and T is temperature. FiglHl illustrates the velocity-force 
characteristics at finite temperature, for different values of T. By increasing T, the effects 
we have shown at T = are increasingly smeared out. In particular, the static threshold 
disappears, as can be seen in the inset of the left panel of FiglHl The bistability regions and 
hysteresis still survive up to small values of T, and the characteristic curve is smoothened 
in the region of dynamical bistability. Interestingly, the area of the hysteresis loop, shown 
in the right panel of FiglHl decreases with T and eventually disappears for sufficiently high 
temperatures. A hysteretic behavior in the intermediate friction region, at T / 0, has 
been reported for long periodic chains Q|. Here we see that we recover essentially the same 
behavior in the case of the dimer. This is because the same mechanisms are effective, as will 
be discussed in the next section. 




(19) 
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FIG. 9: Mobility B =< v > /F, normalized by the asymptotic mobility at large forces Bf = 
(rwy) , as a function of the driven force F. This plot should be compared to that given for 
the Frenkel-Kontorova model by Braun et al. Q| (see their Fig. la). The force is normalized by the 
maximum applied value F max = 0.71kb. The parameters in this case are u = 0.056/c5 2 ,7 = (2/3)«o, 
and a/b = 1/2. 

5. CONCLUDING REMARKS 

Our general aim in the present paper has been to extend previous investigations [3] of 
friction in the simplest non-trivial system, a dimer moving over a periodic substrate, to in- 
clude a driving force and arbitrary temperatures. Our study has shed some light on puzzling 
features in the literature for more complex systems, such as the Frenkel-Kontorova chain, by 
allowing us to obtain simple physical arguments for our system. The model is indeed sim- 
ple: a linear damped oscillator sliding in a sinusoidal periodic potential. Yet, except for the 
limitation that it is restricted to a single spatial dimension, it has the necessary ingredients 
to represent a real dimer or molecule set in a controlled microscopic sliding experiment. For 
example, a molecule sliding along channels of a crystalline well-oriented substrate [lf| should 
exhibit some of the features we have described. Our results, both analytical and numeri- 
cal, confirm the existence of nonlinear friction, resonance effects, bistability, and hysteresis, 
which can be well understood in terms of the resonance of a driven, damped oscillator. Far 
from resonance the sliding friction goes asymptotically to the linear regime with a v ~ 3 term, 
which represents the tail of the resonance response. For intermediate forces, bistability and 
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hysteresis emerge as a consequence of the interplay between linear friction and resonance. 
That our results provide a simple representation of phenomena reported for large systems 
should be clear from FigJHl which shows a comparison between our results on the dimer 
and results from Braun et al. Q for the FVenkel-Kontorova model, which ea„ be regarded 
as the infinite-size generalization of the dimer. Those results are presented to facilitate 
direct comparison with Fig. la of Ref [8] in terms of the mobility B = v/F, normalized 
by its asymptotic value Bf = l/Wy. We notice that the intermediate behavior before the 
asymptotic linear regime is observed in the same fashion as for the dimer. In the latter 
case the mechanism is completely understood as being due to resonance of the dimer; the 
same mechanism could also underlie the nonlinear friction in the Frenkel-Kontorova case 
(see Ref. Jf|). Braun at al. Q attributed such features to the presence of kinks that can be 
observed during the sliding regime of the Frenkel-Kontorova chain. 

In light of the present study, let us make some comments that might have some relevance 
to the ongoing debate in the literature WW about whether S ,id i „ g friction is mai „, y 
electronic or phononic in origin. For this purpose let us identify the background (linear) 
friction in our model with electronic and the resonance friction with phononic sources. Such 
an identification appears natural because, although vibrations of the substrate do not exist 
in our model, resonance friction does arise from the interaction of the dimer vibrations with 
the substrate; on the other hand, the source of the intrinsic friction, denoted in our analysis 
by jv , is independent of interactions with the substrate. We have seen from our analysis that 
the resonance friction is modulated by sin 2 (7ia/b). The commensuration ratio a/b appearing 
in this modulation factor would depend in a realistic 2-d or 3-d environment additionally on 
the relative orientation of the dimer (adlayer) and the substrate. Therefore, while resonance 
friction might dominate the background friction in principle (for sufficiently large corrugation 
amplitude values u), the smallness of the modulation ratio could make it have disparately 
small values relative to background friction. That might be a plausible explanation for the 

nnnn 

disparate results obtained in otherwise similar simulation models (la, 113, 20, 21]. However, 
rather high velocities seem to be necessary for resonance friction to be appreciable. Rough 
estimates we have made suggest that in a number of materials, b ~ 2A, ujq ~ 10 — 100 cm -1 , 
adlayer velocities relative to the substrate necessary for resonance friction to be observable 
would be as high as 30—300 m/s. The velocity region of the quartz microbalance experiments 
corresponds to the low velocity limit discussed at the end of section 13.11 (see EqH2J where 
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the resonance friction becomes linear and the ratio between phononic and electronic friction 
results proportional to (^2) 2 • Let us extend this estimation to the case of a Lennard- 
Jones model for the dimer (adsorbate) interaction, where that ratio becomes proportional 
to (u/e) 2 , e being the well depth of the Lennard- Jones potential. In the Xe over Ag case for 
example e pa 20meV and u can lie in the l-2meV range, therefore the prefactor {u/e) 2 pa 0.01. 
With the above suggested identification of the friction mechanisms — and within the validity 
of analytical assumption used in the low velocity limit — , one might thus expect phononic 
friction to be rarely observable for most materials under typical experimental conditions. 
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